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1 Introduction and the main results 

The problem of finding the irreducible polynomials over finite fields together 
with the related topic of the irreducible factorization of polynomials is one of 
central themes in the theory of finite fields and in the computational algebra, 
and has numerous applications in coding theory, cryptography, computational 
number theory. There are fairly complete surveys of the work in these areas: 
a commehensive account of earlier results can be found in the monograph Q , 
while provides a survey of relatively recent work. The best known algo¬ 
rithm for polynomials irreducible factorization is Berlekamp’s algorithm []|. 
Typically, algorithms for finding irreducible polynomials over a finite field are 
derived from algorithms for irreducible factorization of polynomials provided 
that the irreducible decomposition has a unique factor. A well-known algo¬ 
rithm for determining irreducible polynomials over finite fields is Rabin’s test 
0. There are other similar algorithms (see, for example, i, i). 

In this paper, we propose an algorithm for finding the irreducible poly¬ 
nomials over finite fields, which, in contrast to the aforementioned work, is 
based on the use of the companion matrix of a polynomial, the notion of the 
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multiplicative order of a matrix and the generalized Jordan normal form of a 
linear operator (matrix). 

The companion matrix of a polynomial 

/(t) = ag T 0,it + CL2t^ + + 0,d—2t'^ ^ ^ i 


over a held F is dehned as the matrix 


[/] 


/ 0 1 0 

0 0 1 

0 0 0 

0 0 0 

\-ao -ai -a 2 


0 0 \ 

0 0 

0 0 

0 1 

—ad-2 —ad-i / 


It is not difficult to check that the characteristic polynomial of the companion 
matrix [/] is equal to f{t). 

The multiplicative order (abbreviated to m.o.) of a square matrix A over 
a held F is the minimal positive integer I such that A’’ = E, where E is a unit 
matrix. For the existence of m.o. A it is necessary that detA ^ 0. This condition 
is also a sufficient condition, when F is a hnite held F^, where q = p”, p is 
a prime, n is a positive integer. In this paper only hnite helds F = F^ are 
considered. 

The m.o.[/] of the companion matrix [/] of a polynomial f{t) over Fg with 
a non-zero free term ag is called the order of f{t) and is denoted by ordf{t) 
([1], Dehnition 3.2 and Lemma 8.26). 

Let TO be a positive integer, non-multiple of p. In this (and only in this) 
case there exist positive integers c such that — 1 is divided by to. Any such 
c is divided by the minimal such integer e. This e is called the multiplicative 
order of q{mod.m). 

The main results of this paper are the following theorem and the algorithm, 
based on this theorem. 


Theorem 1 Let F = Fq be a finite field, where q = , p is a prime, n is a 

positive integer. Let f{t) be a unitary polynomial over F with a non-zero free 
term, [f] be its companion matrix and 

d = degf{t), m = ord f{t) = m.o.[f]. (1) 

Then the following assertions are true: 

1) m < q'^ — 1 and m = q’^ — 1 if and only if f{t) is a primitive irreducible 
polynomial; 

2) more generally, f{f) is an irreducible polynomial if and only if p is not a 
divisor of m, d = m.o.qfmod.m) and rk{\fY — E) = d for all positive integers 
I < TO, dividing to. 

Algorithm for finding the irreducible polynomials over a finite 
field 

Suppose that f{t) is a unitary polynomial over a finite field F = Fg of a 
degree d and with a nonzero free term. 
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1. Construct the companion matrix [/] of f{t). 

2. Compute [fY, I = 2,3,... and find m = m.o.[f]- Then m < — 1 and 

m = q‘^ — 1 ii and only if f{t) is a primitive irreducible polynomial. 

3. If TO < — 1 and to is a multiple of p, then the polynomial f{t) is reducible. 

4. Suppose that m < q'^ — 1 and to is not a multiple of p. Then compute 
q\mod.m), I = 2,3,... and find e = m.o.q{mod.m). The polynomial f{t) is 
reducible ii e ^ d. 

5. Suppose that m < q‘^ — 1, m is not a multiple of p and e = d. Then calculate 
r = rk{[fY — E) for divisors h < I 2 < ... of to. If for a divisor ^ < to of 
TO the rank r < d, then /(t) is a reducible polynomial. Otherwise, it is 
irreducible. 

Remark 1. Calculation of degrees of companion matrix [/] can be realized 
in a vertical tape (strip) as follows. Put in the beginning of the strip the matrix 
[/]. At every step, assign a new bottom line, equal to the sum of the last line, 
multiplied by —ad-i, of the penultimate line, multiplied by —ad- 2 , ■■■, of the 
d-th since the end line, multiplied by —qq. Then the square matrix, formed by 
the lines ^ + 1,..., Z + d — 1 of this strip, is equal to [/]*. 

Remark 2. If we know one primitive unitary irreducible over a finite field 
Fg polynomial f(t) of a degree d, then we can find all unitary irreducible over 
Fg polynomials g{t) of a degree d', dividing d, by formula 

g{ty =det{tE Z < to', (Z, to') = 1, (2) 

where d = d's, to' is a divisor of — 1 such that m.o.q{mod.m') = d'. Then 
the order ordg{t) = to'. 

Let us emphasize that by this formula, each unitary irreducible polynomial 
g{t) of degree d' dividing d arises exactly d' times for different values of Z, 
satisfying the condition Z < to', (Z,to') = 1. Thus, we have that the number 
of unitary irreducible over F^ polynomials of a degree d and of an order to, 
which are related by the condition d = m.o.q{mod.in), is equal to (p{m)/d, 
where p{x) is the Euler function. 

Instead of Theorem [T] we prove the following slightly more general result. 

Theorem 2 Let A be a non-singular matrix over a finite field F^ of size dxd, 
f(t) be the characteristic polynomial of A, m = m.o.A. Then the following 
assertions are true. 

1) If f(t) is a reducible polynomial, then 

TO < g'^ - I. (3) 

2) The polynomial f(t) is irreducible if and only if 

(a) TO is not divided by p, 

(b) d = m.o.qfmod.m), 

(c) rk{A^ — E) = d for all I < to, dividing to. 
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Theorem [T] is obtained from Theorem [5] by choosing as matrix A the com¬ 
panion matrix [/] of a unitary polynomial f(t) over F^. Note that the relation 
d = m.o.q{mod.m) implies that m divides q'^ — 1 and that (by definition) f{t) 
is a primitive irreducible polynomial if m = — 1. 


2 Preliminary results 

2.1 Generalized Jordan normal form. 


Our proof of Theorem[2]is based on a theorem about generalized Jordan normal 
form of a matrix (linear operator) (i, i, i)- 

According to this theorem any square matrix A over a field F is conjugated 
to a matrix J = P~^AP, which is a direct sum of generalized Jordan blocks: 

J = Ji©...©Jfe. (4) 


Such a matrix J for any matrix A is determined uniquely up to order of direct 
summands Ji and is called the generalized Jordan normal form of A. 

Any generalized Jordan block Ji is a square block-matrix of a form 



(M w, 0 ... 

0 0 \ 


0 [g^\ W, ... 

0 0 

J^ = 

0 0 [gi\ ... 

0 0 


0 0 0 ... 

[g^] w. 


0 0 0 ... 

0 lgi]J 

where [gi\ is the companion matrix of a unitary irred 


gi{t), which is a divisor of the characteristic polynomial f{t) of A, Wi is a 
square matrix with 1 in the lower left corner and 0 elsewhere. 

For example, if A is the companion matrix [/] of a unitary polynomial f(t) 
and 

f{t) = gi{ty"-9k{ty^ ( 5 ) 


is the irreducible factorization of f{t) with unitary irreducible polynomials 
gi{t), then A has the generalized Jordan normal form Q with Jordan blocks 
Ji, having on the main diagonal si companion matrices of gi(t). 

Since under conjugation the multiplicative order, the rank and the char¬ 
acteristic polynomial of a matrix are preserved, the proof of Theorem [2] is 
reduced to the case, when A is a Jordan matrix 


2.2 Some results on block-matrices. 


We use the following results on block-matrices. 
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The F-algebra of square matrices over the field F of size d x d is denoted 
by Md{F). Let A G Md{F) and d = di + ... + dk, where di,..., dk are positive 
integers. Below we use the following denotations: 

d = (di, ...,dk), Di = di + ... + di,i G fc}. 

Let for arbitrary i,j G k} Aij be the matrix, formed by elements 

of matrix A, which are in intersections of lines Di — di + l,...,Di and rows 
Dj — dj + 1,..., Dj. The matrix A can be considered as a block-matrix with 
blocks Aij, i,j = 1,..., k and d is called the type of this block-matrix. 

Let A,Bg Md{F) be considered as d-block-matrix with blocks Aij and Bij, 
respectively, and C = AB. Then the following property of block multiplication 
holds: the matrix C also can be considered as a d-block-matrix with blocks 

Cij = AiiBij + ... + AikBkj. 

A d-block-matrix A, formed by blocks Aij, is called upper triangular block- 
matrix, if all Aij are zero matrices under i > j. 

In particular, any generalized Jordan block Ji of size di x di, di = Sid'i can 
be considered as an upper triangular d-block-matrix with di = [d'i, ...,d'i). 

Using the property of block multiplication, one can prove that 

1) the product of two upper triangular d-block-matrices A and B with 
blocks Aij and Bij, respectively, is an upper triangular d-block-matrix and its 
diagonal blocks are AuBu; 

2) if A = All © ...© Afcfc, then A™ = Af{ © ...© A™^ for any positive integer 

to ; 

3) if A = All © ... © Akk and each block An has a multiplicative order mi, 
then A has the multiplicative order to = LCM{mi, ...,mk). 


2.3 Two theorems on the multiplicative order of matrices over a finite field. 

Theorem 3 Let J he a generalized Jordan block of a size d x d, which has 
on the main diagonal s companion matrices [g] of a unitary irreducible over 
Fq,q = p” polynomial g{t) of degree d'. Let m' = ordgft) = m.o.[g] and r be 
the minimal integer such that p'’ > d. Then 

m.o.J=p^m'. (6) 

Proof. Considering J as an upper triangular block-matrix of the type d = 
{d',...,d'), we get that is a (usual) upper triangular matrix of a form 
E N, where U is a unit matrix and N is an upper triangular matrix with 
zeros on the main diagonal, having non-zero elements directly after elements of 
the main diagonal. Therefore V is a nilpotent matrix with a nilpotent degree 
d. 

Since the matrices E and N commute, then for any positive integer k we 
have 

{E + = ClE^N^ + ClE^-^N^ + ... + CIe'^N^, 
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where are the binomial coefficients. Since matrices E and N are defined over 
the field F, of characteristic p and the binomial coefficients C], are multiple of 
p under i ^ O^k, when k = p'" for a positive integer r, we get that {E + N)p = 
E + NP\lip'' >d = sd', we have Jp"’"' = = E. 

Suppose that r is the minimal integer such that p^ > d, i. e. 

/ > d, /-I < d. (7) 

Then we obtain the equality (jS]). Theorem |3] is proved. □ 

Using all previous results of this section, we get the following theorem. 

Theorem 4 Let A G Md(Fg) and m = m.o.A. Suppose that A is reduced 
to the generalized Jordan normal form where every generalized Jordan 
block Ji has on the main diagonal Si companion matrices [gi\ of a unitary 
irreducible over polynomials gift) of a degree d) and of a multiplicative 
order m'. Suppose that ri is the minimal integer such that p'"' > di = Sid), if 
Si > 2, and Vi = 0, otherwise. Put r = maxjri,..., r^}. Then 

m = p'"LCMfm'i, (8) 

Note that Theorem 3.11 of Q follows immediately from our Theorem 01 if 
take as matrix A the companion matrix [/] of the polynomial f{t). By using 
the Theorem 3.11 of Q, the inequality m < — 1 of Theorem [1] of this 

paper can be obtained from the inequality p'^LCM{m{, — 1. In 

subsection 3.2 below this inequality is proved in a more general context of our 
Theorem 01 


2.4 On binomials t‘> — t. 

For proving the basic results of this paper we also need the following theorem. 
Theorem 5 Any binomial 

bq,d(t) =t‘^'‘ -t 

is the product of all unitary irreducible polynomials over the field F^ of a degree 
d', dividing d, each with the multiplicity one. 

Proof. In its field of decomposition F {bq^d{t)), the binomial bq^d{t) has exactly 
q‘^ different roots, because bf ^{t) = —1. The set of roots of the binomial bq^d{t) 
is closed under operations of subtraction and division on non-zero elements, 
consequently, this set coincide with the field F {hq^dft)). 

Let g{t) be an irreducible over F^ polynomial of a degree d'. Then the field 
Fq[t[l{g{t)) consists of elements, consequently, it is the field of decomposi¬ 
tion for the binomial bq^d'it). The following four assertions are equivalent: 

(i) the polynomial gft) divides the binomial bq^d(t); 

(ii) the field F {bq^d'{t)) is isomorphic to a subfield of the field F {hq^d{t))] 

(hi) — 1 divides — 1; 

(iv) d' divides d. 

Therefore the assertions of Theorem 0 are true. □ 
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3 Proofs of assertions of Section 1 

3.1 Proof of assertion 2 of Theorem [5] 

Firstly, starting the proof of Theorem [21 we suppose that the characteristic 
polynomial f{t) of a matrix A G Md{Fq) with g = p” is irreducible, and prove 
that then m = m.o.A is not a multiple of p and d = m.o.qijnod.m). 

If /(t) G Fg[l] is an arbitrary irreducible polynomial of degree d and (/(t)) 
is the prime ideal, generated by f{t), then Fq[t]/{f{t)) is a finite field with q'^ 
elements and its multiplicative group is a cyclic group of order q‘^ — 1. 

Recall that according to Cayley-Hamilton theorem, the matrix A is a root 
of f{t). Since by supposition f{t) is an irreducible polynomial, the extension 
Fg(A) of the field F^ by A can be considered as a matrix model of the field 
Fq[t]/{f{t)) in the ring of matrices Md{Fq). Therefore, according to Lagrange 
theorem, m = m.o.A is a divisor of q'^ — 1. Hence 

(i) m is not divided by p and 

(ii) d is an exponent for q[mod.m). 

Since p does not divide m, there exists e = m.o.q{modm), and since d is 
an exponent for q(mod.m), then e divides d. 

On the other hand, the relation e = m.o.q{mod.m) implies that m divides 
q^ — If TO = m.o.A, we get that = E. Thus, H is a root of the 

binomial Since, according to Cayley-Hamilton theorem, the irreducible 

(characteristic) polynomial /(t) also has A as a root, then f{t) is decomposed 
in the field Fq(A) = F {bq^d{t)) on the linear factors and, consequently, f{t) 
divides hq^eit)- Then, according to Theorem jS] d divides e. 

Thus, we have proved that to is not divided by p and d = e = n.o.q{mod.m). 
To finish the proof of the direct implication of the second assertion of Theorem 
121 left to prove that rk{A^ — E) = d, \i I < m and I divides to. 

We have that is a root of the binomial bq^d{t). Then according to The¬ 
orem [5] it is a root of an irreducible polynomial g{t), dividing bq^d(t), with a 
degree d' = degg{t), dividing d. 

Suppose that d = d's. Then the characteristic polynomial of A* is equal to 
g{tY and the generalized Jordan normal form of A^ is a direct sum of s copies 
of companion matrices [p]. Therefore rk{A^ — E) = d, \i I < m. 

Now we prove the inverse implication of the second assertion of Theorem 
121 i.e. the sufficiency the conditions (a), (b), (c) of Theorem[21for irreducibility 
of the characteristic polynomial f{t) of matrix A. According to the subsection 
2.1, it is sufficient to check it for Jordan matrices 01 

From the condition (a) and Theorems [3l and |4l we have that each Jordan 
block Ji is the companion matrix of an irreducible polynomial gi(t). 

From the condition (c) we get that all gi{i) have the sane multiplicative 
order to' and, consequently, a same degree 

d' = m.o.q(mod.m'). (9) 


Hence according to 01 


d = d'k 


(10) 
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and according to Theorem H] 


m = m.o.J = m'. (11) 

Then, using the condition (b) and the relations (Pl. lfTUl) . (ITTl) . we get that 
fc = 1. The second assertion of Theorem [5] is proved. 


3.2 Proof of assertion 1 of Theorem [5] 

In this subsection we check the inequality ([3]). Actually, it is sufficient to check 
this inequality for generalized Jordan matrices J = Ji © ... © Jfc. 

If the unitary polynomial f{t) over is reducible over this field, then 
either f{t) = g{tY, where g{t) is a unitary irreducible polynomial over 
and the integer s > 2 , or f{t) = fi{t)f 2 {t), where fi{t) and / 2 (<)are different 
unitary (therefore non-constant) polynomials. Correspondingly, there are two 
following alternative possibilities for the Jordan normal form J of A: 

(i) J consists of a unique generalized Jordan block, which has on the main 
diagonal s > 2 companion matrices [ 5 ] of an irreducible polynomial g(t); 

(ii) J is a direct sum of fc > 2 generalized Jordan blocks. 

We check the inequality m for these two cases separately. 

Consider the case (i). Then J is a matrix of the size d x d, where d = 
sd', d' = deg g{t). 

By using Theorem [3] and the inequalities 

1 ) p’’ < d of the condition (O, 

2) m' < q'^ —1, which is true, because we have proved in item 2 that the multi¬ 
plicative order m' of the matrix [( 7 ] with the irreducible over Fg characteristic 
polynomial g{t) of order d' divides q'^ — !> ( or because d' = m.o.q{mod.m')), 

3) P < g = p", 

we obtain that the inequality ([3]) follows from (is equivalent to) the inequality 

pd's < + {p'^'y-^ +...+ p‘‘' + 1 , ( 12 ) 

Under s > 3 this inequality is true, because px < p^, if a; = 1,2, and 
p <px < p“, if X > 3. 

Under s = 2 the inequality m, evidently, also is true, excepting the cases 

(a) d' = 1 and p is an arbitrary prime; 

(b) d' = 2 and p = 2 or 3; 

(c) d' = 3 and p = 2. 

Thus, the inequality m is proved, excepting these three cases. Below we 
check that the inequality Q holds in these cases as well. 

(a) We have d' = 1, s = 2, d = sd' = 2 and from (O we get that r = 1. 
Therefore the inequality JS]) is reduced to the inequality p < p + 1. 

(b) d' = s = 2, d = 4 and p is 2 or 3. Then the condition (I7|) implies r = 2. 
Consequently, the inequality m take a form p^ < p^ + 1 . 
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(c) d' = 3, s = 2,d = 6,p = 2. The condition d?]) gives r = 3 and the 
inequality m take a form + 1. 

Thus, the inequality ([3]) for matrices J with a reducible characteristic poly¬ 
nomials f{t) is proved in the case (i). 

Consider the case (ii). Let J be a direct sum of k Jordan blocks Ji of sizes 
di X di and of multiplicative orders rrii, i = 1,k. For any Jordan block we 
have proved above that m.o.Ji = rrii < — 1. Using the block form for J of 

type d= (di,...,dfc), we get that 






Consequently, 


T.J < — 1). 




If s > 2, we have 

- l)...{q‘^*‘ - 1) < q^ - 1, d = di + ... + 4 

Thus, the first assertion of theorem [3] and, consequently, all this theorem are 
completely proved. 


3.3 Proofs of assertions of Remarks 1 and 2. 

The assertion of Remark 1 follows immediately by the definition of multipli¬ 
cation of matrices. So, it remains to prove the assertions of Remark 2. 

Let f{t) be a unitary primitive irreducible over the finite field F, polyno¬ 
mial of a degree d. Then [ff, I = 1,..., — 1 together with the corresponding 

zero matrix are all roots of the binomial bq^dit). On the other hand, the bino¬ 
mial bq^d{t) is equal to the product of all irreducible over the field F, unitary 
irreducible polynomials g(t) of a degree d', dividing d, with a multiplicity one. 
Therefore each such g{t) has d' roots 

[f]^\ I <m', {l,m') = 1 

and, naturally, is uniquely determined by this set of roots. Here 

m' = ordg{t) = m.o.[g], d' = m.o.q{mob.m'). 

Then each such matrix has a characteristic polynomial g{tY and its general¬ 
ized Jordan normal form is equal to the direct sum of s copies of companion 
matrices [g]. This explicates Remark 2. 
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